Recovery of a Medieval Brucella melitensis genome using shotgun metagenomics

ABSTRACT Shotgun metagenomics provides a powerful assumption-free approach to the recovery of pathogen genomes from contemporary and historical material. We sequenced the metagenome of a calcified nodule from the skeleton of a 14th-century middle-aged male excavated from the medieval Sardinian settlement of Geridu. We obtained 6.5-fold coverage of a Brucella melitensis genome. Sequence reads from this genome showed signatures typical of ancient or aged DNA. Despite the relatively low coverage, we were able to use information from single-nucleotide polymorphisms to place the medieval pathogen genome within a clade of B. melitensis strains that included the well-studied Ether strain and two other recent Italian isolates. We confirmed this placement using information from deletions and IS711 insertions. We conclude that metagenomics stands ready to document past and present infections, shedding light on the emergence, evolution, and spread of microbial pathogens

Identification of novel imidazo[1,2-a]pyridine inhibitors targeting M. tuberculosis QcrB.

Mycobacterium tuberculosis is a major human pathogen and the causative agent for the pulmonary disease, tuberculosis (TB). Current treatment programs to combat TB are under threat due to the emergence of multi-drug and extensively-drug resistant TB. Through the use of high throughput whole cell screening of an extensive compound library a number of imidazo[1,2-a]pyridine (IP) compounds were obtained as potent lead molecules active against M. tuberculosis and Mycobacterium bovis BCG. The IP inhibitors (1–4) demonstrated minimum inhibitory concentrations (MICs) in the range of 0.03 to 5 µM against a panel of M. tuberculosis strains. M. bovis BCG spontaneous resistant mutants were generated against IP 1, 3, and 4 at 5× MIC and subsequent whole genome sequencing identified a single nucleotide polymorphism 937ACC>937GCC (T313A) in qcrB, which encodes the b subunit of the electron transport ubiquinol cytochrome C reductase. This mutation also conferred cross-resistance against IP 1, 3 and 4 demonstrating a common target. Gene dosage experiments confirmed M. bovis BCG QcrB as the target where over-expression in M. bovis BCG led to an increase in MIC from 0.5 to >8 µM for IP 3. An acute murine model of TB infection established bacteriostatic activity of the IP series, which await further detailed characterization

Genomic diversity of Escherichia coli isolates from backyard chickens and guinea fowl in the Gambia

Chickens and guinea fowl are commonly reared in Gambian homes as affordable sources of protein. Using standard microbiological techniques, we obtained 68 caecal isolates of Escherichia coli from 10 chickens and 9 guinea fowl in rural Gambia. After Illumina whole-genome sequencing, 28 sequence types were detected in the isolates (4 of them novel), of which ST155 was the most common (22/68, 32 %). These strains span four of the eight main phylogroups of E. coli, with phylogroups B1 and A being most prevalent. Nearly a third of the isolates harboured at least one antimicrobial resistance gene, while most of the ST155 isolates (14/22, 64 %) encoded resistance to ≥3 classes of clinically relevant antibiotics, as well as putative virulence factors, suggesting pathogenic potential in humans. Furthermore, hierarchical clustering revealed that several Gambian poultry strains were closely related to isolates from humans. Although the ST155 lineage is common in poultry from Africa and South America, the Gambian ST155 isolates belong to a unique cgMLST cluster comprising closely related (38-39 alleles differences) isolates from poultry and livestock from sub-Saharan Africa - suggesting that strains can be exchanged between poultry and livestock in this setting. Continued surveillance of E. coli and other potential pathogens in rural backyard poultry from sub-Saharan Africa is warranted

Remarkable genomic diversity among Escherichia isolates recovered from healthy chickens

The genus Escherichia has been extensively studied and it is known to encompass a range of commensal and pathogenic bacteria that primarily inhabit the gastrointestinal tracts of warm-blooded vertebrates. However, the presence of E. coli as a model organism and potential pathogen has diverted attention away from commensal strains and other species in the genus. To investigate the diversity of Escherichia in healthy chickens, we collected fecal samples from antibiotic-free Lohmann Brown layer hens and determined the genome sequences of 100 isolates, 81 of which were indistinguishable at the HC0 level of the Hierarchical Clustering of Core Genome Multi-Locus Sequence Typing scheme. Despite initial selection on CHROMagar Orientation medium, which is considered selective for E. coli, in silico phylotyping and core genome single nucleotide polymorphism analysis revealed the presence of at least one representative of all major clades of Escherichia, except for E. albertii, Shigella, and E. coli phylogroup B2 and cryptic clade I. The most frequent phylogenomic groups were E. coli phylogroups A and B1 and E. ruysiae (clades III and IV). We compiled a collection of reference strains isolated from avian sources (predominantly chicken), representing every Escherichia phylogroup and species, and used it to confirm the phylogeny and diversity of our isolates. Overall, the isolates carried low numbers of the virulence and antibiotic resistance genes typically seen in avian pathogenic E. coli. Notably, the clades not recovered are ones that have been most strongly associated with virulence by other studies

Evolution of Mycobacterium tuberculosis complex lineages and their role in an emerging threat of multidrug resistant tuberculosis in Bamako, Mali

In recent years Bamako has been faced with an emerging threat from multidrug resistant TB (MDR-TB). Whole genome sequence analysis was performed on a subset of 76 isolates from a total of 208 isolates recovered from tuberculosis patients in Bamako, Mali between 2006 and 2012. Among the 76 patients, 61(80.3%) new cases and 15(19.7%) retreatment cases, 12 (16%) were infected by MDR-TB. The dominant lineage was the Euro-American lineage, Lineage 4. Within Lineage 4, the Cameroon genotype was the most prevalent genotype (n=20, 26%), followed by the Ghana genotype (n=16, 21%). A sub-clade of the Cameroon genotype, which emerged ~22 years ago was likely to be involved in community transmission. A sub-clade of the Ghana genotype that arose approximately 30 years ago was an important cause of MDR-TB in Bamako. The Ghana genotype isolates appeared more likely to be MDR than other genotypes after controlling for treatment history. We identifed a clade of four related Beijing isolates that included one MDR-TB isolate. It is a major concern to fnd the Cameroon and Ghana genotypes involved in community transmission and MDR-TB respectively. The presence of the Beijing genotype in Bamako remains worrying, given its high transmissibility and virulence

Lateral gene transfer between prokaryotes and multicellular eukaryotes: ongoing and significant?

The expansion of genome sequencing projects has produced accumulating evidence for lateral transfer of genes between prokaryotic and eukaryotic genomes. However, it remains controversial whether these genes are of functional importance in their recipient host. Nikoh and Nakabachi, in a recent paper in BMC Biology, take a first step and show that two genes of bacterial origin are highly expressed in the pea aphid Acyrthosiphon pisum. Active gene expression of transferred genes is supported by three other recent studies. Future studies should reveal whether functional proteins are produced and whether and how these are targeted to the appropriate compartment. We argue that the transfer of genes between host and symbiont may occasionally be of great evolutionary importance, particularly in the evolution of the symbiotic interaction itself

Evolution of virulent genotypes and an emerging threat of multidrug resistant tuberculosis in Bamako, Mali

Background: Bamako, Mali, has a well-structured hierarchyfor tuberculosis (TB) case management. However, in recent yearsBamako has been faced with an emerging threat from multidrugresistant TB (MDR-TB). Here we present insights into the genomicepidemiology of TB and the evolutionary mechanisms driving theemergence of MDR-TB in Bamako.Methods & Materials:Isolates recovered from tuberculosispatients from local reference centers and the University Teach-ing Hospital at Point G, in Bamako, Mali between 2006 and 2012(n = 208), were tested for antimicrobial susceptibility at the MRCUnit The Gambia. A subset of 76 isolates were analysed using wholegenome sequencing. A time dated phylogenetic tree was recon-structed using BEAST. Lineage and resistance conferring mutationswere inferred using PhyResSe.Results:Patients included 21 females and 55 males agedbetween 3 to 78 years, among whom 12(16%) were infectedby MDR-TB. Most patients 61(80%) were new cases and among15 retreatment cases 9(60%) were MDR-TB. The phylogeny wasreconstructed from 8508 variant core genome sites. The dom-inant lineage was the Euro-American super lineage, lineage 4.Within lineage 4, the Cameroon genotype was the most preva-lent genotype (n = 20, 26%) followed by the Ghana genotype (n = 16,21%). Cameroon genotype isolates diverged from a common recentancestor∼161 years ago to form three clusters, one of whichemerged∼22 years ago and is likely to be involved in on-goingtransmission. Seven Ghana genotype isolates were MDR-TB repre-senting over half all MDR-TB in this dataset (7/12). Ghana genotypeisolates were more likely to cause MDR-TB than other genotypesafter controlling for treatment status (OR = 5.6, p-value = 0.043).The MDR-TB Ghana genotype isolates formed a clade that divergedapproximately 30 years ago, in which thekatGSer315Thr mutationwas conserved. Other Euro-American genotypes included the sixLAM, two H37Rv-like and one Uganda. Four patients were infectedwith closely related Beijing strains and five patients were infectedwith non-MDRMycobacterium africanum2

A new approach to in silico SNP detection and some new SNPs in the Bacillus anthracis genome

<p>Abstract</p> <p>Background</p> <p><it>Bacillus anthracis </it>is one of the most monomorphic pathogens known. Identification of polymorphisms in its genome is essential for taxonomic classification, for determination of recent evolutionary changes, and for evaluation of pathogenic potency.</p> <p>Findings</p> <p>In this work three strains of the <it>Bacillus anthracis </it>genome are compared and previously unpublished single nucleotide polymorphisms (SNPs) are revealed. Moreover, it is shown that, despite the highly monomorphic nature of <it>Bacillus anthracis</it>, the SNPs are (1) abundant in the genome and (2) distributed relatively uniformly across the sequence.</p> <p>Conclusions</p> <p>The findings support the proposition that SNPs, together with indels and variable number tandem repeats (VNTRs), can be used effectively not only for the differentiation of perfect strain data, but also for the comparison of moderately incomplete, noisy and, in some cases, unknown <it>Bacillus anthracis </it>strains. In the case when the data is of still lower quality, a new DNA sequence fingerprinting approach based on recently introduced markers, based on combinatorial-analytic concepts and called cyclic difference sets, can be used.</p

GlyGly-CTERM and Rhombosortase: A C-Terminal Protein Processing Signal in a Many-to-One Pairing with a Rhomboid Family Intramembrane Serine Protease

The rhomboid family of serine proteases occurs in all domains of life. Its members contain at least six hydrophobic membrane-spanning helices, with an active site serine located deep within the hydrophobic interior of the plasma membrane. The model member GlpG from Escherichia coli is heavily studied through engineered mutant forms, varied model substrates, and multiple X-ray crystal studies, yet its relationship to endogenous substrates is not well understood. Here we describe an apparent membrane anchoring C-terminal homology domain that appears in numerous genera including Shewanella, Vibrio, Acinetobacter, and Ralstonia, but excluding Escherichia and Haemophilus. Individual genomes encode up to thirteen members, usually homologous to each other only in this C-terminal region. The domain's tripartite architecture consists of motif, transmembrane helix, and cluster of basic residues at the protein C-terminus, as also seen with the LPXTG recognition sequence for sortase A and the PEP-CTERM recognition sequence for exosortase. Partial Phylogenetic Profiling identifies a distinctive rhomboid-like protease subfamily almost perfectly co-distributed with this recognition sequence. This protease subfamily and its putative target domain are hereby renamed rhombosortase and GlyGly-CTERM, respectively. The protease and target are encoded by consecutive genes in most genomes with just a single target, but far apart otherwise. The signature motif of the Rhombo-CTERM domain, often SGGS, only partially resembles known cleavage sites of rhomboid protease family model substrates. Some protein families that have several members with C-terminal GlyGly-CTERM domains also have additional members with LPXTG or PEP-CTERM domains instead, suggesting there may be common themes to the post-translational processing of these proteins by three different membrane protein superfamilies

Sequence Diversities of Serine-Aspartate Repeat Genes among Staphylococcus aureus Isolates from Different Hosts Presumably by Horizontal Gene Transfer

BACKGROUND: Horizontal gene transfer (HGT) is recognized as one of the major forces for bacterial genome evolution. Many clinically important bacteria may acquire virulence factors and antibiotic resistance through HGT. The comparative genomic analysis has become an important tool for identifying HGT in emerging pathogens. In this study, the Serine-Aspartate Repeat (Sdr) family has been compared among different sources of Staphylococcus aureus (S. aureus) to discover sequence diversities within their genomes. METHODOLOGY/PRINCIPAL FINDINGS: Four sdr genes were analyzed for 21 different S. aureus strains and 218 mastitis-associated S. aureus isolates from Canada. Comparative genomic analyses revealed that S. aureus strains from bovine mastitis (RF122 and mastitis isolates in this study), ovine mastitis (ED133), pig (ST398), chicken (ED98), and human methicillin-resistant S. aureus (MRSA) (TCH130, MRSA252, Mu3, Mu50, N315, 04-02981, JH1 and JH9) were highly associated with one another, presumably due to HGT. In addition, several types of insertion and deletion were found in sdr genes of many isolates. A new insertion sequence was found in mastitis isolates, which was presumably responsible for the HGT of sdrC gene among different strains. Moreover, the sdr genes could be used to type S. aureus. Regional difference of sdr genes distribution was also indicated among the tested S. aureus isolates. Finally, certain associations were found between sdr genes and subclinical or clinical mastitis isolates. CONCLUSIONS: Certain sdr gene sequences were shared in S. aureus strains and isolates from different species presumably due to HGT. Our results also suggest that the distributional assay of virulence factors should detect the full sequences or full functional regions of these factors. The traditional assay using short conserved regions may not be accurate or credible. These findings have important implications with regard to animal husbandry practices that may inadvertently enhance the contact of human and animal bacterial pathogens